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Abstract 

We propose a novel RG method to specify the location of the IR fixed point in lattice gauge theories and apply it to the 
517(3) gauge theories with Nf fundamental fermions. It is based on the scaling behavior of the propagator through the 
RG analysis with a finite IR cut-off, which we cannot remove in the conformal field theories in sharp contrast with the 
confining theories. The method also enables us to estimate the anomalous mass dimension in the continuum limit at the 
IR fixed point. We perform the program for Nf = 16,12,8 and Nf = 7 and indeed identify the location of the IR fixed 
points in all cases. 


Scale invariance, or more precisely conformal invari¬ 
ance has become a fundamental concept in understanding 
the universal aspects of the nature from the Planck scale 
to the Hubble scale. They appear not only in critical phe¬ 
nomena of condensed matter physics, but also in quantum 
gravity, high energy particle phenomenology, and all the 
way up to cosmology [l[. Many conformal field theories, 
however, are strongly coupled, and much remains unsolved 
in their theoretical understanding. In particular, when re¬ 
alized by gauge theories, the constructive approaches to 
the conformal fixed points are still rudimentary [2]. The 
aim of this article is to clarify some important aspects of 
these constructive approaches and offer one simple crite¬ 
rion on conformal invariance. 

Obviously, the central question is to locate the IR fixed 
point within a given class of theories. In this article we pro¬ 
pose a novel and simple RG method to specify the location 
of the IR fixed point in lattice gauge theories by studying 
the scaling behavior of the propagator. We will apply the 
technique to the SU( 3) gauge theories with Nf fundamen¬ 
tal fermions (within the conformal window), and estimate 
the anomalous mass dimension. We perform this program 
for Nf = 16,12,8 and Nf = 7, and indeed identify the 
location of the IR fixed points in all cases. 

We constructively define gauge theories on Euclidean 
plane R 4 as the continuum limit of lattice gauge theories 
on the Euclidean lattice of the size N x = N y = N z = N 


and Nt = rN (r being an aspect ratio, which is fixed as 
r = 4 throughout the article), taking the limit of the lattice 
space a —> 0 and N —» cso, with L = N a and L t = N t a 
fixed. When L and/or L t are finite, the system is bounded 
by an IR cutoff Air ~ 1 /L. We impose an anti-periodic 
boundary condition in the time direction for fermion fields 
and periodic boundary conditions otherwise. In conformal 
field theories the IR cutoff is an indispensable ingredient 
because there is no other natural scale to compare, which 
will be further elucidated in this article. 

Our general argument that follows can be applied to 
any gauge theories with fermions in arbitrary (vector-like) 
representations, but to be specific, we focus on 51/(3) 
gauge theories with Nf fundamental Dirac fermions. For 
the lattice regularization of the action, we employ the Wil¬ 
son quark action and the RG improved gauge action Q 
(also known as the Iwasaki gauge action in the literature). 

Given the regularized action, the theory is defined by 
two parameters; the bare coupling constant go and the bare 
degenerate quark mass mo at ultraviolet (UV) cutoff. We 
also use, instead of go and mo, (3 = 6/<?q and the hopping 
parameter K = 1/2 (mod + 4). 

As for observables, together with the plaquette and the 
Polyakov loop in each space-time direction, we measure the 
quark mass m q defined through Ward-Takahashi identities 


(OIV 4 A 4 IPS) 
q 2(0|P|PS) 


(1) 
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where P is the pseudo-scalar density and A 4 the fourth 
component of the local axial vector current, with renor¬ 
malization constants being suppressed. The quark mass 
Triq defined in this way only depends on /3 and K up to 
order 1/7V corrections. 

One of the most important observables we will study 
is the t dependence of the propagator of the local meson 
operator in the H channel: 

G H (t) = , (2) 

X 

where the summation is over all the spatial lattice points. 
In this paper, we mostly focus on the pseudo-scalar (PS) 
channel H = PS , and the subscript H is suppressed here¬ 
after. 

In order to investigate the large t behavior of a propa¬ 
gator, we define the effective mass m(t) through 

cosh (rn{t){t - N t /2)) _ G(t) 
cosh(m(t)(f + 1 — N t /2)) G(t + 1) 

When boundary effects can be neglected, it reduces to 

m(t)=1 "Gfriy (4) 


transformation induced by the change of the UV renor¬ 
malization scale 3 ' = fi/s, followed by a space-time scale 
change by a factor 1 /s (see e.g. 0), relates the propaga¬ 
tor with different parameters as 

fN'\ 3-27 

G{t;g,m q ,N,ii) = f— J G{t\g' ,m' q ,N' ,/x). (6) 

Here N' = N/s and t' = t/s. The relation between g' and 
g and m' q and m q are determined by the beta function B 
and the mass anomalous dimensions 7 . 

The UV renormalization scale 3 in lattice theories is 
set by the inverse lattice spacing a -1 . In Eq. © the pa¬ 
rameter g, is common on both sides and therefore g, may 
be suppressed in the notation, but we should keep in mind 
that the RG equation ([ 6 ]) implicitly assumes the contin¬ 
uum limit V, N' —y 00 to avoid the effect of the UV cutoff. 

Let us first discuss the case in which we are at the fixed 
point, i.e. g' = g = g* and m' q = m q = 0 so that B = 0 and 
7 = 7 *. In this case, the propagator may have simplified 
notation as 

G(r,N)=G(t,N). (7) 

with r = t/N t . The variable t takes 0, 1 , 2, • • • , N t — 1 so 
that 0 < t < 1. In terms of r, the RG relation eq. © 
reduces to 


In the case of exponential-type decay the effective mass 
approaches a constant value in the large t regime, which 
we call a plateau. 

Before non-perturbative discussion, let us first recall 
the perturbative result. Within the two-loop perturbation 
theory, the RG beta function for the SU (3) gauge coupling 
constant is given as 


B(g) = - 


(33 — 2N f ) 3 
487r 2 


(102-f N f ) 


(167T 2 ) 2 


9 b + 0(g 7 ) . 


(5) 


The fixed point B(g*) = 0 exists for 8.05 < Nf < 
16.5 within the two-loop approximation!!]]. When Nf = 
16, the IR fixed point is located at /? = 11.48. Since this 
coupling constant is small, we may trust the perturbative 
computations for Nf = 16. We will compare the non- 
perturbative calculation with this value. 

When Nf decreases, g* increases, at least in the per¬ 
turbation theory, and therefore non-perturbative effects 
become important. The smallest Nf where the fixed point 
exists is denoted as Nf and the range of flavors Nf < Nf < 
16 is called the “conformal window”. The lower bound 
of the conformal window can only be determined non- 
perturbatively. Our earlier studies [5!]@] strongly suggest 
the conjecture that the conformal window is 7 < Nf < 16. 
However, the conjecture is based on indirect logics. In this 
article, we will present more direct evidence supporting the 
conjecture. 

Let us study the RG properties of the propagator in the 
vicinity of the fixed point. The RG equation for the RG 


/ at' \ 3-27 * 

G(r;V)=(—j G(t;N). ( 8 ) 

Strictly speaking, this equation is satisfied in the limit 
N , N' ->■ 00 . 

To state our proposal concretely, we define the scaled 
effective mass m(£; N) as 

“ |l ’ l,| = l,h cpu») ' (9) 

In the continuum limit N —> 00 Eq. © reduces to the 
form 

m(r, N) = — d T In G(r, N) (10) 

The crucial observation, which will be the core of our pro¬ 
posal is that, combining Eqs.© and (flOll . the scaled effec¬ 
tive mass does not depend on N as a function of r: 

m(r, N) = m(r, N ) (11) 

at the fixed point. Therefore, the agreement of the scaled 
effective mass as a function of N and r are stringent tests 
of the fixed point. 

Suppose that we are away from the fixed point (i.e 
g =£ g* while m q = 0) in contrast. The scaled effective 
mass in the vicinity of the fixed point would instead show 
the following behavior in the leading term 

m( t, g, N) = m(r, g, N') + B(g) In 5 g m(r, 3 , N') . 

( 12 ) 
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Effective mass: Nf=16; beta=11.0, K=0.129 


Effective mass: Nf=16; beta=10.0, K=0.1294 




Figure 1: Scaled effective mass plots for Nf = 1 6 at. 8 = 11.0 and 10.0: three sets of symbols are N = 16 (red square), N = 12 (green circle), 
N = 8 (blue triangle). 


Effective mass: Nf=16; beta=10.5, K=0.1292 



t/N t 


Effective mass: Nf=16; beta=10.5, K=0.1292 



t/N, 


Figure 2: Scaled effective mass plots for Nf = 16 at ft = 10.5: the left panel is an enlarged one of the right panel; three sets of symbols are 
N = 16 (red square), N = 12 (green circle), N = 8 (blue triangle). 


We can rewrite Eq. in a more suggestive form 
m(r, g, N) - m(r, g, N') 


m = 




(13) 


which may be used to compute B(g) in the vicinity of the 
fixed points once we numerically compute m(r, g , N). This 
formula is satisfied in the limit N, N 1 — > oo. 

Our strategy is as follows. With given Nf and /3, we 
tune the quark mass (defined through Ward-Takahashi 
identity) to be zero. Then we numerically compute the 
meson propagator on the lattice. For each choice of the 
lattice size N, we plot the effective mass defined by Eq. 
© in terms of the scaled time r. As we explained, gener- 
ically, the scaled effective mass do not coincide with each 
other as a function of t at a given value of /3 but different 
values of N. However, if we find the fixed point value /3, 
the plots for different N must coincide with each other. 


In this article, we perform numerical simulations on the 
three lattices with size 8 3 x 32,12 3 x 48 and 16 3 x 64 with 
the aspect ratio of r = 4. Note that although the explicit 
form of m(r, N) depends on r, the relation Eq. (fill) itself 
is valid for any r. 

Let us show in Fig. 1 the scaled effective mass plots 
in two cases of many such examples: we take Nf = 16 
at /3 = 11.0 and /3 = 10.0. The asymptotic behaviors of 
three sets of data points and the lines connecting them on 
N = 8,12 and 16 lattices, do not coincide with each other. 
We may conclude that these values of /3 do not correspond 
to the fixed point. On the other hand, as we will see in Fig 
2, if we take /3 = 10.5, then the three plots and the lines 
do coincide within the standard deviation. Based on the 
RG relations, we claim that this is the value of the gauge 
coupling constant at the fixed point. 

We perform this program for Nf = 7, 8,12,16 on lat- 
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Table 1: Simulation parameters for the cases we identify the IR fixed 
points: the first column N is the lattice size, the second column -/V tra j 
is the number of trajectories, the third ”acc” is the acceptance ratio 
for the global Metropolis test, ”plaq” is the value of the plaquette 
and m q is the quark mass defined by eq. Q. 



TV/ = 

= 16, /T = 

= 10.5, K = o.: 

1292 

TV 

TV-traj 

acc 

plaq 

rriq 

16 

2000 

0.59(1) 

0.922 55(1) 

-0.0063(1) 

12 

4000 

0.77(1) 

0.922 55(1) 

-0.0053(1) 

08 

4000 

0.89(1) 

0.922 57(1) 

0.0003(5) 


Nf 

= 12,/3 

= 3.0, K = 0.1405 

TV 

TVtraj 

acc 

plaq 

m q 

16 

3000 

0 . 68 ( 1 ) 

0.74416(2) 

- 0 . 002 ( 1 ) 

12 

3000 

0.84(1) 

0.74415(1) 

- 0 . 002 ( 1 ) 

08 

4000 

0.94(1) 

0.74419(2) 

0.004(1) 


£ 

II 

00 

= 2.4, K = 0.147 

TV 

TVtraj 

acc 

plaq 

m q 

16 

4000 

0.72(1) 

0.676 20(1) 

-0.007(1) 

12 

4000 

0.84(1) 

0.676 20(1) 

-0.006(3) 

08 

3000 

0.93(1) 

0.676 22(2) 

-0.0005(5) 


Nf 

= 7,/3 = 

= 2.3, I< = 0.14877 

TV 

TVtraj 

acc 

plaq 

m q 

16 

4000 

0.72(1) 

0.659 31(1) 

-0.0017(2) 

12 

4000 

0.85(1) 

0.659 31(1) 

-0.0005(3) 

08 

5000 

0.94(1) 

0.65941(3) 

0.0047(6) 


tices with size 8 3 x 32,12 3 x 48 and 16 3 x 64. By narrowing 
down the region where the scaled effective mass m(r, TV) 
becomes close for different TV, we identify the point /?* 
where they agree with each other within one standard de¬ 
viation. Note that the scaled effective mass as a function 
of r in general does depend on TV, so the coincidence of 
the three curves only at a particular common value of /3 is 
dramatic. We eventually find that this occurs for all the 
cases we study. 

The algorithms we employ are the blocked HMC algo¬ 
rithm Q| in the case Nf = 2 N and the RHMC algorithm 
0] for Nf = 1 in the case Nf = 2 N + 1. We show the 
parameters of simulations for the cases we identify the IR 
fixed points in Table I. The masses of quarks are expressed 
in units of a -1 . We choose the run-parameters in such a 
way that the acceptance of the global Metropolis test is 
about 60% ~ 90%. The statistics are 1,000 MD trajecto¬ 
ries for thermalization and 1000 ~ 4000 MD trajectories 
for the measurement. We estimate the errors by the jack¬ 
knife method with a bin size corresponding to 100 HMC 
trajectories. 

Before examining our numerical results, we have a few 
comments in order. 

Firstly, for gauge configuration generation we have to 
be very careful to choose the lowest energy state, not 
quasi-stable states in the conformal region. As stressed 
in Ref. [ 6 ], there are quasi-stable states, which persist for 
long time for a HMC algorithm. 


Secondly, in contrast with the confining phase, when 
the system is either in the deconfining phase at high tem¬ 
perature or in the conformal region, it is not hard to per¬ 
form simulations at zero quark mass. It is even possi¬ 
ble to calculate across the zero quark mass from posi¬ 
tive to negative mass without any trouble. This is be¬ 
cause in the deconfining phase the density of eigenvalues of 
the massless Dirac-Wilson operator decreases toward zero 
('modification [l2j| [l3| of the Banks and Cacher relation^.) 
We used this fact to identify the first order chiral phase 
transition for Nf = 3 and 6 , which we call “on the K c 
method” |Tol |. and to find the fact that for Nf > 7 there 
is no confining phase at the massless quark in the strong 
coupling limit. 

Thirdly, because the chiral symmetry is explicitly bro¬ 
ken in the Wilson action, we have to tune the hopping 
parameter K. The quark mass does have 1/N correction. 
The mass at N = 12 differ from that at N = 16 with order 
of 0(0.001), while that at TV = 8 with order of 0(0.005) in 
our case. We estimate the effect on the meson propagator 
of this difference. The effect by the difference of 0(0.001) 
is one of order smaller than the statistical errors and that 
of 0(0.005) is order of a half of one standard deviation. In 
total, we estimate the smallness of the difference is enough 
for the accuracy we take in this article. 

Now, let us show the results, starting with the Nf = 16 
case. As mentioned earlier, within the two-loop perturba¬ 
tion, the IR fixed point is /T* = 11.48, which is RG scheme 
independent. On the other hand, the coupling constants 
in different RG schemes are related to each other by a 
constant as /3i = fa + c \2 in the one-loop approximation. 
For example fill] . the lattice coupling constants /Trg and 
/Tone—plaquette for one-plaquette action are related to that 
in the continuum theory /T^g (in the modified minimal 
subtraction scheme) as 


/Trg = /3ms - 0.3 


and 


/Tone—plaquette — /3 ms A 3.1. 


It is well-known that the convergence of the perturbation 
by the g 0 ne-piaquette is poor in general. The contribution 
of higher order terms will be large. On the other hand, 
the lattice coupling constant /3rg is close to 13^0 [14| and 
therefore we may expect that the higher-order contribution 
is not so large and the location of the fixed point is close 
to 11.2 in /Trg from the two-loop estimate and the above 
relation. 

In order to find the fixed point from our proposal,we 
make several trials including those at /3 = 10.0,10.5,11.0, 
and 11.5. We find the three sets of data and the lines 
connecting them are apparently different from each other 
at /3 = 11.0 (Fig. 1; left panel) and they approach closer by 
decreasing /3 as 11.0,10.5. On the other hand, at /T = 10.0 
(Fig. 1; right panel) they are apart each other again but 
they approach closer by increasing /J as 10.0,10.5. This 
suggests that there is an IR fixed point between /T = 10.0 
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Effective mass: Nf=12; beta=3.0, K=0.1405 


Effective mass: Nf=08; beta=2.4, K=0.147 




t/N, t/N, 


Figure 3: Scaled effective mass plots for Nf = 12 at 0 = 3.0 and Nf = 8 for /3 = 2.4; three sets of symbols are N = 16 (red square), .V — 12 
(green circle) and N = 8 (blue triangle). 


Effective mass: Nf=07; beta=2.3, K=0.14877 



t/N, 


Figure 4: Scaled effective mass plots for Nf = 7 at 8 = 2.3: three 
sets of symbols are N = 16 (red square), N = 12 (green circle) and 
N = 8 (blue triangle). 


and 11.0. We indeed find, as shown in Fig. 2, that the three 
sets of the scaled effective mass plots are almost degenerate 
at /3 = 10.5 and K = 0.1292. We see that three lines almost 
overlap for r > 0.1. Only in the small r region (r < 0.1) we 
see the differences. We interpret the difference for r < 0.1 
is due to the fact that N is not large enough to remove 
the effect of the UV cutoff /r = a -1 . 

The fact that our method identifies the location of the 
IR fixed point at a value expected from the perturbation 
theory, together with the fact that three lines almost over¬ 
lap, strengthens our confidence in the validity of our ap¬ 
proach. 

We make similar process for Nf = 12,8 and 7 as the 
Nf = 16 case. In Figs. 3 and 4 are shown the results. 
The qualitative feature of our results are the same. If we 
choose a very particular /3 for each Nf , the data and three 


Propagator: Nf=8; beta=2.4, k=147 



t/N, 


Figure 5: Propagators for Nf = 8 at /3 = 2.4: three sets of symbols 
are N = 16 (red square), N = 12 (green circle) and N = 8 (blue 
triangle). 


lines almost overlap for r > 0.1, as shown in the Figures. 
In the small t, region (t/N t < 0.1) we find the differences. 
Since they are similar to the case of Nf = 16, we do not 
present them here. 

Finally we identify the IR fixed points at /3* = 10.5±0.5 
for Nf = 16; 3.0 ± 0.2 for N f = 12; 2.4 ± 0.1 for N f = 8; 
and 2.3 ± 0.05 for N f = 7. 

On the other hand, in the Nf = 6 case, there is a chiral 
phase transition point at finite (3 when N is finite [Ic|- If we 
would perform a program similar to the above (by fixing 
/3 and increasing the lattice size N ), then at some N the 
system would end up with the confining phase rather than 
the chiral symmetric phase (to which the conformal fixed 
points belong). Thus the IR behavior would be completely 
different. It cannot be a conformal field theory. 

Thus our results at the finite lattice size (up to 16 3 x64) 
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are consistent with that the conformal window is 7 < Nf < 
16. However we do not exclude the possibility of the “walk¬ 
ing scenario” that the RG beta function is anomalously 
small near the edge of the conformal window (e.g. Nf = 7 
or 8 ), and for a larger N an undiscovered chiral phase tran¬ 
sition point happens to appear at some value of /3 and the 
chiral phase transition eventually occurs in the infinite N 
limit. 

Eq. (]i:ip may be used to verify whether the sign of the 
beta function changes at the fixed points identified, calcu¬ 
lating the beta function just above and below the critical 
coupling constant. The formula contains the differences 
of the scaled effective masses in the denominator and nu¬ 
merator. The differences are of order 0(1)% of the scaled 
effective masses and the errors of our data are also of the 
same order. Therefore it is hard to estimate numerically 
the beta function just above and below the critical cou¬ 
pling constant with our present data. 

The simulation at larger lattices will clarify which sce¬ 
nario is realized in the continuum limit, verifying first 
whether the fixed point does stay at the identified fixed 
point and second whether the sign change occurs. 

We note that it seems possible to extract the mass 
anomalous dimension 7 * using the scaling of propagators 
Eq.(© in the continuum limit N — > 00 . In Fig [5] we show 
the results for the propagator on the three lattice sizes in 
the Nf = 8 case. The data are depicted on a logarithmic 
scale. We see the data roughly scale at r > 0.1. Similar 
results are obtained also for Nf = 16,12 and Nf = 7. 

We try to estimate 7 * using the scaling law in Eqj 8 ] 
For each Nf there are three sets for N and N': (8,16), 
( 8 , 12), (12,16). In each set we estimate 3.0 — 27 * and 
statistical errors at each corresponding r, treating data at 
each r as statistically independent. The statistical errors 
are small at small r and become large at large r, so we 
quote the results at r = 0.125. The results are in order of 
(8,16), ( 8 , 12 ), (12,16): 

N f = 16:- 0.0309(24), -0.0516(24), -0.0017(64). 

N f = 12 :0.0337(35), 0.0086(55), 0.0692(105). 

N f = 8 :0.15148(51),0.1047(77), 0.1289(125). 

Nf = 7 :0.1767(55), 0.1652(83),0.1929(163). 

In the case Nf = 16, the perturbative computation gives 
7 * = 0.026. The order of magnitude for Nf = 16 is rea¬ 
sonable given the large statistical error. The order of es¬ 
timated anomalous mass dimensions for Nf = 12 is ap¬ 
parently smaller than those reported in the literature (See 
for example E3). The reason may be N and N' are not 
large enough to obtain anomalous mass dimensions in the 
continuum limit. 

To gain an intuition about the systematic error, we 
can perform the same scaling analysis for the free massless 
Wilson fermion. It shows the “anomalous dimension” of 
—0.06, —0.07, —0.04 at r = 0.125 (for the same lattice size 
ratio above) while the scaling of effective masses is less 
affected by the finite size effects. In general, the effective 


masses scale better than the amplitude even with smaller 
N. To find the complete resolution of the issue we have to 
simulate with lager N. We leave it for the future study. 

Let us finally discuss the physical meaning of the shape 
of the effective mass plot at the fixed point. In order to 
investigate the dynamics of the theory, we proposed a de¬ 
tailed analysis of temporal propagators which we call the 
“local-analysis” of propagators [ 6 |. We parametrize the 
propagator G(t) as 


G{t) = c(t) 


exp(— m(t) t) 


(14) 


It is possible to determine c(to), m(to), a(to) locally, using 
three point data G(to),G(to + 1), G(t 0 + 2 ). This is not a 
fit, but a parametrization. In addition to this parametriza- 
tion, we also perform a fit using data at five points and 
check the fit and the parametrization are in good agree¬ 
ment, which implies that the parametrization represents 
the characteristics of the propagator reasonably well. 

The results of the local parameters are similar to those 
obtained in our previous article @ except for some quanti¬ 
tative differences. Our current data actually suggests that 
the data in our previous workQ are slightly off the IR 
fixed points. Nevertheless almost all what were discussed 
there can be applied also here. To avoid repetition we do 
not plot m(t), but plot a(t) for Nf = 16,12,8, and 7 in 
Figs. 02 and [7] since the a(t) mostly represents the char¬ 
acteristics. We observe an interesting change of the form 
of a(t) from Nf = 16 to Nf = 7. For Nf = 16 there is 
a shoulder around a ~ 1.5 at 8 < t < 14. On the other 
hand for Nf = 8 and 7, there is a plateau around a ~ 0.7 
at t > 14. The Nf = 12 case shows a transition between 
them. 

In addition to the local analysis of propagators we ex¬ 
amine the Polyakov loops in spatial directions to investi¬ 
gate the vacuum structure. The results are also similar to 
those in Ref. [f|. As a typical example we show the scat¬ 
tered plot of the Polyakov loops forTV/ = 16 in Fig (8] 

From these results we may conclude that the vacuum 
of the Nf = 16 is the Z( 3) twisted one modified by non- 
perturbative effects and the meson is an almost free fermion 
state in the twisted vacuum. The t dependence of a{t) is 
very similar to that of free fermions in the Z{ 3) twisted 
vacuum. See Ref. [ 6 ] for the details. 

The plateaus at 14 < t < 24 in the Nf = 8 and 7 cases, 
taking values a(t) ~ 0.7 ± 0.1 is well described by an un¬ 
particle meson model and 7 * ~ 1.3 ± 0.1 from the formula 
2 — a = 7 * 0. We need the data with high statistics for a 
more precise value of 7 *. 

We would like to stress that a conformal field theory is 
completely different from QCD in the point that there is 
no dimensional parameter such as Aqcd- In QCD if Na is 
large enough compared with Aqcd, boundary effects can 
be neglected and it can be assumed the limit TV = 00 is 
taken. However the boundary effects are essential even 
at any large lattice TV in the conformal field theories be¬ 
cause there is no other natural scale to compare. Note 
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Nf=16, Beta=10.5, K=0.1292 Nf=12, Beta=3.0, K=0.1405 




Figure 6: Local exponent a(t): Nf = 16 in the left panel; Nf = 12 in the right panel 


Nf=8, Beta=2.4, K=0.147 Nf=7, Beta=2.3, K=0.14877 




Figure 7: Local exponent a(t): Nf = 8 in the left panel; Nf = 7 in the right panel 


that the propagators Eqs. (7) and (11) are functions of 
the scaled time r which takes value 0.0 < r < 1.0. Clearly 
the function depends on the boundary condition as well as 
the aspect ratio even if we take N —> oo limit. Of course, 
to be clear, this does not mean that the local physics of 
the conformal field theory depends on the boundary con¬ 
ditions we use. We note the zero momentum propagator 
in our definition m may not be a local variable because 
we have summed over spatial coordinates before taking the 
continuum limit. 

In the near future we would like to perform the program 
with larger lattice sizes and more statistics to derive the 
anomalous mass dimension using Eqj51 and the relation of 
the eigenvalue density of the Dirac-Wilson operator and 
the anomalous index 121 Tdj ljj . It would be intriguing to 
compare them with the value from the unparticle meson 
model. 

The calculations were performed with HA-PACS com¬ 
puter at CCS, University of Tsukuba and SR16000 at 


KEK. We would like to thank members of CCS and KEK 
for their strong support for this work. 
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